Discrete skew self adjoint canonical systems 
and the isotropic Heisenberg magnet model 

M.A. KAASHOEK and A.L. SAKHNOVICH 
Abstract 

A discrete analog of a skew selfadjoint canonical (Zakharov-Shabat 
or AKNS) system with a pseudo-exponential potential is introduced. 
For the corresponding Weyl function the direct and inverse problem 
are solved explicitly in terms of three parameter matrices. As an 
application explicit solutions are obtained for the discrete integrable 
nonlinear equation corresponding to the isotropic Heisenberg magnet 
model. State space techniques from mathematical system theory play 
an important role in the proofs. 

Introduction 

In this paper we shall treat a discrete analog of the well-known skew selfad- 
joint canonical (Dirac type, Zakharov-Shabat or AKNS) system: 



dY 

-iJ——(x,z) = zY(x,z) + V(x)Y(x,z), J 
ax 



Im 

-/„ 



, x > 0. (0.1) 



Here z is a spectral variable, Y and V are 2m x 2m matrix functions on the 
half line, and V is skew selfadjoint, that is, V(x)* = JV(x) with V(x)* being 
the matrix adjoint of V(x). To obtain the discrete analog of (|U.1|) let U be 
the unique solution of the initial value problem 

^-(x) = -iU(x)JV(x), x>0, U(0) = I 2m . (0.2) 
ax 
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Since JV(x) is selfadjoint, we get from (|U.2|) that U(x) is unitary for each 
x > 0. Now put S(x) = U{x)JU{x)* and W(x,z) = U{x)Y(x,z). Then 

^-(x, z) = izS(x)W(x, z), S(x) = S(x)* = S(xy\ x>0. (0.3) 
ax 

It is now immediate that 

W n+1 (X)-W n (X) = -jS n W n (X), S n = S* n = S~\ n = 0, 1,... (0.4) 

is a natural discrete analog of 1)0.1)1 . This discrete analog of the continuous 
pseudo-canonical system is very important. In fact, when m = 1, then the 
system ()0.4)) turns out to be an auxiliary system for the nonlinear isotropic 
Heisenberg magnet (IHM) model [3B] (see also the detailed discussion after 
Theorem 10.51 below and the historical remarks in [121 )■ Motivated by the 
IHM model we shall use the term spin sequence to denote any sequence of 
N x N matrices {S n } satisfying 

S n = S* n = S-\ n = 0,1,2,... (0.5) 

As for the skew selfadjoint continuous case [22] (see also [HI]), one can 
associate with ()0.4j) an m x m matrix function (p(X), meromorphic on < 
— 5 < 0, such that 

im]w n {xyw n {x) 

n=0 

where W n (X), n > 0, is the fundamental solution of (j0.4)) . i.e., the 2m x 2m 
matrix solution W n (X) of ()0.4j) normalized by the condition Wo(X) = I 2m . 
One refers to (p as the Weyl function of ()0.4j) . When the Weyl function is 
rational and strictly proper we shall recover the system ()0.4)) explicitly from 
its Weyl function. For this purpose we need to introduce spin sequences that 
are the discrete analogs of the pseudo-exponential potentials from JTHJ HZ] 
(see also the references therein). 

The spin sequences from this special class are defined in terms of three 
parameter matrices with the following properties. First fix an integer N > 0, 
and consider anJVxiV matrix a with det a ^ 0, an N x N matrix S such 



cp(X) 



< oo, 



(0.6) 
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that E = Eq, and an N x 2m matrix A . These matrices should satisfy the 
following matrix identity 

aE - E a* = z'A Aq. (0.7) 

Given these three matrices a, So, and Ao we define for n = 1,2,... the 
N x 2m matrix A n and the N x N matrix E n via recursion: 

A n+ i = A n + ia' l k n J, 

E n+1 = E n + o ; - 1 E ri (a*)- 1 + Q ; - 1 A n ^A;(Q ; *)- 1 . (0.8) 

Next assume that the matrices E n , n = 0,1,2,..., are non-singular. Then 
we say that the sequence of matrices {S n } defined by 

S n = J + A^E^-A^E^Am-i, n = 0,1,2,..., (0.9) 

is the spin sequence determined by the parameter matrices a, Eo and Ao- 
Notice that this requires the invertibility of the matrices E n . 

For spin sequences defined in this way our first theorem presents an 
formula for the fundamental solution W n (X) of (|0.4|) . 

Theorem 0.1 Let a (det a ^ 0), E (S = Eq) and A satisfy <\U.1\) , and 
assume that det E n ^ for < n < M, where E n is given by For 
< n < M — 1 let S n be the matrices determined by a, E and A via ()^.^|) 
and (123). Then S n = S* = S' 1 for < n < M - 1, and for0<n<M the 
fundamental solution W^A) of the discrete system (\0.4l can be represented 
in the form 

/ i \ n 

W n (X) = W a>A (n,\)(l 2m - -JJ W aA (0, A)" 1 , (0.10) 

where W a ^{n, A) is defined by 

W aA (n, A) = I 2m + zA^E^A/at - a)- l K. (0.11) 

When E > 0, there exist simple conditions on a and A to guarantee 
that det E n ^ 0. First, if So > 0, then without loss of generality we can 
assume that E = In- Indeed, it is easy to see that the sequence of matrices 
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{S n } defined by (|U.9|) and (|U.8|) does not change if we substitute a, S and 

_i i _i 

Ao by S 2 q;Sq, In and E 2 Ao- So let us assume that So = In- Next, we 
partition A into two N x m blocks 9i and #2 as follows: A = [#1 9 2 \- This 
together with S = In allows us to rewrite ()0.7|) in the form 

a -a* = i(9x9l + 9 2 e;). (0.12) 

Furthermore, in this case A n is given by 

A n = [(I N + ia- 1 ) n 9 1 {lN-iu- l ) n 9 2 }. (0.13) 

Finally, we shall assume that the pair {a, 9i} is full range which means that 

= span{a%C m | k = 0, 1, 2, . . . , N - 1}. 

The following proposition shows that under these conditions automatically 
det a^O and det H n ^ for n — 0, 1, 2, 

Proposition 0.2 Zet a fc o square matrix of order N , and 9\ and 6*2 be 

N x m matrices satisfying tyO.lty . Assume that the pair {a, 9\} is full range. 
Then all the eigenvalues of a are in the open upper half plane C + , and for 
n = 1,2,... the matrices S n defined by §0.8ty . with S = In an d A n given by 
(J7T73J), are positive definite and satisfy the identity 

aS n - S n a* = iA n A* . (0.14) 

Definition 0.3 A triple of matrices a, 9\ and 92, with a square of order N 
and 9\ and ^2 of size N xm, is called admissible if the pairs {a, 9i} and {a, 
9 2 } are full range and the identity $0. lty holds. 

We denote by the acronym FG (finitely generated) the class of spin sequences 
{S n } determined by the matrices a, S = I N and A = [6>i 9 2 ], where a, 
9\ and 9 2 form an admissible triple. In this case we also say that these spin 
sequences are determined by the corresponding admissible triples. The next 
two theorems present the solutions of the direct and inverse problem in terms 
of the Weyl function. 
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Theorem 0.4 Assume that the spin sequence {S n } n >o of the discrete pseudo- 
canonical system f\U.4l belongs to the class FG and is determined by the ad- 
missible triple a, 9\ and 9 2 . Then the system (HO) has a unique Weyl function 



(p, which satisfies (jfl. 6|) on the half plane 3 A < — |, a finite number of poles 
excluded, and this function is given by the formula 

i f (\)=t6* 1 (\I N -{3)- 1 9 2 , (0.15) 

where 

f3 = a-i9 2 9* 2 . (0.16) 

Notice that the function tp in 1)0.15)1 is a strictly proper m x m rational 
matrix function. Conversely, if ip is a strictly proper m x m rational matrix 
function, then it admits a representation of the form 

ip{\) = i$l{\I n - 1 )- 1 $ 2 , (0.17) 

where 7 is a square matrix and $i,$2 are matrices of size n x m. We refer 
to the right hand side of ([0.17)1 as a minimal realization of <p if among all 
possible representations ()0.17j) of p the order n of the matrix 7 is as small as 
possible. This terminology is taken from mathematical system theory. We 
can now state the solution of the inverse problem. 



Theorem 0.5 Let ip be a strictly proper rational m x m matrix function, 
given by the minimal realization < \U. 1T[ . There is a unique positive definite 
n x n matrix solution X of the algebraic Riccati equation 

7X-X7* = i(X$ l d* 1 X -d 2 d* 2 ). (0.18) 

Using X define matrices 81,62, and a = (3 + iQ^Q^ 

9 1 =xH u e 2 = X-*# 2 , (3 = X-^X^. (0.19) 

Then a, 61, and 9 2 form an admissible triple, and the given matrix function 
ip is the Weyl function of a system (\0.4l of which the spin sequence {S n } G 
FG and is uniquely determined by the admissible triple a, 9\, and 9 2 . 
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Next, we describe connections with the nonlinear IHM equation. For this 
purpose consider the zero curvature representation of the IHM model: 

j t G n (t, A) = F n+1 (t, X)G n (t, A) - G n (t, X)F n (t, A), (0.20) 

where 

G n (t,X) = I 2 -^S n (t), F n (t,X) = ^- + ^-, (0.21) 
A A — i A + 1 

V?(t) := (1 + ~3 r -i(t) ■ 'Sr(t))-\h ± S r (t))(I 2 ± S r -!(t)). (0.22) 

Here the vectors S r — [S^ S% Sf] belong to M 3 , R is real axis, the dot ■ 
denotes the scalar product in R 3 , and the correspondence between the spin 
matrix S r and the spin vector S r is given by the equality 



Sr 



C3 Cl A C2 



(0.23) 



In other words the IHM equation 



dSn 25 r „/\( + J- 1 _ ), (0.24) 

' v Vl+ S n ■ S 



dt 



where /\ stands for the vector product in IR 3 , is equivalent to the compatibility 
condition ()0.20|) of the systems 

W n+1 (t,X) = G n (t,X)W n (t,X), 

^-W n (t,X) = F n (t,X)W n (t,X) (n>0). (0.25) 
at 

In fl0.24|) it is required that 1 1 S n \ \ = 1 . Now one can see easily that the 
representation ()0.23|) . where S n ■ S n — 1, is equivalent to the equalities 
(10. 5|) with S n 7^ ±I 2 . Thus the first system in (I0.25|) coincides with system 
f!0.4|) . where m = 2, and S n ^ ±I 2 . We use these connections to obtain 
explicit solutions of the IHM model. 

The literature on continuous canonical systems is very rich, especially 
for the selfadjoint case; see, for instance, the books (TUJ E21 UH 1213 EH]- 
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Selfadjoint continuous canonical systems with pseudo-exponential potentials 
have been introduced in ^3] , and for this class of potentials various direct and 
inverse problems have been solved; see JJJ and the references therein. The 
subclass of strictly pseudo-exponential potentials has been treated in [3IUE]- 
Interesting recent results on the spectral theory of selfadjoint discrete systems 
and various useful references on this subject can be found in El EH I2H1 
137] . Mainly Jacobi matrices (or block Jacobi matrices as in that are 
related to Toda chain problems have been studied. For the skew selfadjoint 
discrete case some references can be found in |TT) IT2"] ITS] . 

Theorem lU.ll is the discrete analog of Theorem 1.2 in The right hand 
side of (jU.llj) can be viewed as the transfer function of a linear input-output 
system (see [6 ). Transfer functions of the special form given by (jO.llj) were 
introduced in [33], and also used for the representation of the fundamental 
solutions of continuous canonical systems [33J |33] . In Theorem IU.1I we are 
closer to jTB] (see also j3Ul|32|), where the dependence on the parameter x 
differs from the one in (33J|3J)]. The condition on an admissible triple a, d%, 
82 that the pairs {a, 81} and {a, 82} are full range pairs is specific for the 
discrete case. Nevertheless Theorems II) .4[ 10.51 and parts of their proofs are 
analogous to results and proofs in Section 2 of [TB] . 

This paper consists of four sections not counting this introduction. Since 
elements from mathematical system theory play an important role in this 
paper, we present in the first section the necessary preliminaries from that 
area. In the second section we prove Theorem 10 . 1 1 and present some auxiliary 
results that will be used in the application to the IHM model. Theorems 10.41 
10.51 and Proposition 10.21 are proved in Section 3. In Section 4 we construct 
solutions of the IHM equation (|0.24j) . describe the evolution of the Weyl 
function and consider a simple example. 

1 Preliminaries from mathematical system 
theory 

The material from the state space theory of rational matrix functions, that 
is used in this paper, has its roots in the Kalman theory of input-output 
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systems |2T], and can be found in books, see, e.g., [121 IE]- In general the 
rational matrix functions appearing in this paper are proper, that is, analytic 
at infinity, and they are square, of size m x m, say. Such a function F can 
be represented in the form 

F(X) = D + C(XI N -A)- 1 B, (1.1) 

where A is a square matrix (of which the order iV may be much larger than 
m), the matrices B and C are of sizes N x m and m x N, respectively, and 
D = F(oo). In this paper D is often a zero matrix, and in that case F is 
called strictly proper. The representation (jl.lj) is called a realization or a 
transfer matrix representation of F, and the number N = ord(A), is called 
the state space dimension of the realization. Here ord(A) denotes the order 
of the matrix A. 

Realizations of a fixed F are not unique. The realization (jl.lj) is said 
to be minimal if its state space dimension iV is minimal among all possible 
realizations of F. The state space dimension of a minimal realization of F 
is called the McMillan degree of F and is denoted by degF. Notice that 
degF = corresponds to the case when ord(A) = 0, and this occurs if and 
only if F(X) = D. The realization (jl.lj) of F is minimal if and only if 

N-l N-l 

span [J Im A k B = C N , f] Ker CA k = {0}, N = oid{A). (1.2) 

k=0 k=0 

If for a pair of matrices {^4, B} the first equality in (jl.2j) holds, then {^4, B} 
is called controllable or a /«// range pair. If the second equality in (jl.2j) is 
fulfilled, then {C, A} is said to be observable or a zero kernel pair. If a pair 
{A, B} is full range, and K is an m x N matrix, where iV is the order of A 
and m is the number of columns for B, then the pair {A — BK, B} is also 
full range. An analagous result holds true for zero kernel pairs. 

Minimal realizations are unique up to a basis transformation, that is, 
if (jl.lj) is a minimal realization of F and F(A) = D + C(XIn — A)~ 1 B is 
a second minimal realization of F, then there exists an invertible matrix S 
such that 

A = SAS-\ B = SB, C = CS-\ 
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In this case S is called a state space similarity. 

Finally if in (jl.ljl we have D = I m , then F(X) is invertible whenever A is 
not an eigenvalue of A — BC and in that case 

F(X)- 1 = I m -C{\I N - A x y l B, A X =A-BC. (1.3) 

2 The fundamental solution 

In this section we prove Theorem IU.1I and present some results that will be 
used in Section 0] (a result on the invertibility of matrices S n , in particular). 

Proof of Theorem 10. 1L First we shall show that equalities ()0.7|) and ()U.8|) 

yield the identity (j(J.14j) for all n > 0. The statement is proved by induction. 
Indeed, for n = it is true by assumption. Suppose (j(J.14j) is true for n = r. 
Then using the expression for £ r +i from (|0.8J) and identity (|0.14|) for n = r 

we get 

o;S r _|_i — E r _|_ia* = 
= iA r A* r + ia^KKWY 1 + A^A^aT 1 ~ a" x A r JA r *. (2.1) 

The first relation in ()0.8|) and formula ()2.1j) yield ()0.14|) for r = n + 1 and 
thus for all n > 0. 

The next equality will be crucial for our proof. Namely we shall show 
that for < n < M — 1 we have 

W aA (n + I, A)(/ 2m - i J) = (I 2m - jS n )W a , A (n, A). (2.2) 

By (|0.11|) formula ()2.2|) is equivalent to the formula 

j(S n -J) = (hm-^A^iXlN-a^An- 

-A^-^XIn - aY^+iihm - jJ). (2.3) 
Using the Taylor expansion of (XI n — at infinity one shows that ()2.3j) 
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is in its turn equivalent to the set of equalities: 

S n — J = A*S n 1 A„, - A* +1 E n+1 A n+ i, (2.4) 

= A;£;> p A n - ^ n A;S; 1 a p - 1 A n (p > 0). (2.5) 

Equality ()2.4j) is equivalent to ()0.9j) . Taking into account the first relation in 
()0.8j) we have <xA n+ i — ik n+ iJ = aA n + a _1 A n . Thus the equalities in (|2.5jl 
can be rewritten in the form i^ n a p_2 A n = 0, where 

K n = A^-^a 2 + I N ) - KK 1 ® 2 + iSnKPn 1 *- (2-6) 

Therefore, if we prove that K n = 0, then equalities ()2.5j) will be proved, and 
so formula ()2.2|) will be proved too. Substitute (jl).9J) into (|2.6p . and again 
use the first relation in ()0.8|) to obtain 

K n = A* n+1 J:-l 1 (a 2 + I N )-A* n J:- 1 a 2 + iJA* n i:- 1 a + 

+^A;S; 1 A n A;S; 1 a - zA; +1 S-^(A n + ia~ l A n J) A*^ 1 a. {2.1) 

Now we use (|U.14|) to obtain iA^A^S" 1 = a — S„a*S~ 1 and substitute this 
relation into ()2.7|) . After easy transformations it follows that 

K n = A; +1 S^ 1 (a- 1 S n (a*)- 1 + S n + a- 1 A n JA;(a*)- 1 )a*S; 1 a + 

+iJA;S; 1 a-AyE;'a (2.8) 

In view of the second relation in (jU.8j) the first term on the right hand side 
of (|2.8|) equals A* +1 a*S~ 1 a and we have 

K n = (A* +1 + iJA* n (a*)- 1 - A* n )a*^- l a. (2.9) 

By the first relation in (|U.8|) the equality K n = is now immediate, i.e., (|2.2|) 
is true. 

Notice that equality (j0.10|) is valid for n = 0. Suppose that it is valid for 
n = r. Then ijOl) and (HHU1) yield 

Wr+i(A) = (J 2m - ^ r )W aj A(r, A)(/ 2m - ^j) V q ,a(0, A)- 1 . (2.10) 
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By and (gHHD the validity of (jlTTUJ) for n = r + 1 easily follows, i.e., 
(ID.lUp is proved by induction. 

Consider now the matrices S n given by ()0.9|) . It is easy to see that 
S n = S*. Notice also that in view of (j0.14|) we have 

W aA (r,X)W a , A (r,\y = I 2m (r > 0), (2.11) 

where A stands for complex conjugate for A. From ()2.2)1 and ()2.11|) it follows 
that (I 2m - i\~ 1 S n )(l2m + iX^Sn) = A _2 (A 2 + l)h m - Thus the equality 
S* = S' 1 holds, which finishes the proof of the theorem. □ 

The case when ±i (jL a(a) (a means spectrum) is important for the study 
of the IHM model. Assume this condition is fulfilled, and put 

R n = (I N - ia- l y n Y, n {lN + ^aT 1 )-™, (2-12) 

Q n = (I N + ia~ l )- n Y, n {I N ~ i{a*r l )- n . (2.13) 

The following proposition will be useful for formulating the conditions of 
invertibility of S n in a somewhat different form then those in Proposition l0.2l 
(see Corollary 12. 21 below) . 

Proposition 2.1 Let the matrices a (det a ^ 0), S = Sq, and A satisfy 
!\U. ?D , and let the matrices S n be given by Ifi ^ c(a), then the sequence 

of matrices {R n } is well defined and non- decreasing. If —i ^ cr(a), then the 
sequence of matrices {Q n } is well defined and non-increasing. 

Proof. To prove that the sequence {R n } is non-decreasing it will suffice to 
show that 

S„+i - (I N - ia- 1 )^^ + i(a*y 1 ) > 0. (2.14) 
For this purpose notice that 

S„+i - (I N - ia _1 )S n (/ A r + i(a*) _1 ) = 

= S n+ i - S n - a' l T, n (a*Y 1 - ia~ l (all n - S n a*)(a*) _1 . 

Hence, in view of ()0.8|) and ()0.14J) we get 

S„+i - (I N ~ ia _1 )S n (J7v + i(a*y r ) = 

= a~\k n JK + KK)(°?)~ X - (2-15) 
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Since J + J 2m > 0, the inequality (|2.14|) is immediate from (|2.15|) . 
Similarly, from (|U.8|) and f|(J. 14|) we get 

E n+1 - (I N + ia -1 )£ n ( I N -i{a*)~ x ) = 

= a-\A n JA* n - A n A;)(a*)- 1 < 0,(2.16) 

and so the sequence of matrices {Q n } is non-increasing. □ 

According to Proposition ^. 11 when i £ a(a) and £ > 0, we have R n > 0. 

Corollary 2.2 Let the conditions of Proposition [Ql hold, and assume that 
i ^ a(a) and So > 0. Then we get E n > for all n > 0. 

Partition the matrices W a ^{r, A) and A r into two m-column blocks each: 

W a , A (r, A) = [{W a , A (r, X)) 1 {W a , A (r, A))J , A r = [(A r ) x (A r ) 2 ]. 
The next lemma will be used in Section 4. 

Lemma 2.3 Let the matrices a (0, ±z ^ cr(ct)), S (S = £q) and Aq satisfy 
l\U. 7| ) ; and /ei t/ie matrices E n 6e gwen 6?/ (Ifl.&j) . Then for n > the following 
relations hold: 

( W Q , A (n, z)) a = (W a ,A(« + 1, x 

x (/ m + 2(W a , A (n, ijJjA^;^ 2 + / JV )- 1 (A n ) 1 ), (2.17) 

( W a , A (n, -i)) 2 = (W aiA (n + 1, i)) 2 x 

x (/ m - 2(W a , A (n, -z))*A;E^ V + J iV )- 1 (A n ) 2 ) . (2.18) 

Proof. From the proof of Theorem 10. II we know that K n = 0, where K n is 
given by (J2.6)) . In particular, we get 

K n ar\c? + I N )-\A n ) x = Q, K n a-\a 2 + In)- 1 (A n ) 2 = 0. (2.19) 

To prove ()2. 17|) notice that (A n+ i)i = a~ 1 (a+U]sr)(A n )i and rewrite the first 
equality in (|2.19|) as 

A; +1 E^ x (q + iI N )-\A n+1 ) 1 - A^ia - z/ Af )- 1 (A n ) 1 + 

+i{hm + S n )A*Jl- l {a 2 + / JV )- 1 (A n ) 1 = 0. (2.20) 
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Put A = — i in (|2.2|) and take into account (|2.11j) to derive 

I 2m + S n = 2(W aA (n + 1, -fy^W^n, i))* r (2.21) 
In view of definition (|D~TT|) of W aA , equality (|2~T7j) follows from (l2~2TH) and 

Putting A = i in ()2.2|) we get 

/ 2m - S n = 2(W aA (n + l 7 i)) 2 (W a , A (n, -i))* 2 . (2.22) 

Analogously to the proof of ()2.17|) we derive ()2.18|) from ()2.22|) and the second 
equality in (ETTHj) . □ 

Remark 2.4 According to and (\2. H| the rank of the matrices l2m±S n 

is less than or equal to m. Together with the formula $U.5\ this implies that 
under the conditions of Lemma we have S n = U n JU*, where U n are 
unitary matrices and J is defined in $0. ij) . 

3 Weyl functions: direct and inverse prob- 
lems 

In this section we prove Theorems IU.4I and IU.51 and Proposition IU.21 At the 
end of the section a lemma on the case i G" a(a) is treated too. 

Proof of Proposition 10.21 Suppose / is an eigenvector of a, that is, 
af = cf, f 7^ 0. Then formula ()0.12|) yields the equality 

i(c - c )/7 = f*(e 1 e* 1 + e 2 e*)f > o. (3.1) 

So c G C + . Moreover, if c = c, then according to (|3.1|) we have 6^f = 6^f = 0, 
and therefore af — a*f. It follows that 

f*e 1 = 0, f*(a-cI N ) = (f^O). (3.2) 

As {a, 6{\ is a full range pair, so the pair {a — cljy, 9{\ is full range, which 
contradicts ()3.2j) . This implies that c G C+, i.e., a(a) C C + . 
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Recall that identity (j(J.14j) was deduced in the proof of Theorem ID. II 
Taking into account that a (a) C C+, identity (|(J.14|) yields 

1 f°° 

S n = — / (a- XI^KKia* ~ XlN)- l d\. (3.3) 

Notice now that the pair {a, + ia~ l ) n 6i\ is full range and use (|U.13|) . 
(Q to obtain S n > for all n > 0. □ 

Remark 3.1 In the same way as in the proof of Proposition \0.S\ above the 
inclusion o~(a) C C + follows from the weaker condition that the pair {a, A } 
is full range. However, the example N — 1, a — i, Q\ — ; #2^2 — 2, which 
yields S n = for n > 0, shows that we have to require that the pair {a, #1} 
is full range in order to get E n > 0. The full range condition on the pair {a, 
A } is not enough for this conclusion. 

Recall now Definition 10 .3l of the admissible triple. Proposition l0.2l implies. 
in particular, that det a 7^ for the admissible triple and the spin sequences 
{S n } determined by it are well defined for all n > 0. In other words the class 
FG is well defined. 

Proof of Theorem IOI Let W a , A (n, A) be given by flUTEJ) . Write W Q)A (0, A) 
as 

a(A) b(X) 



W a , A (0,X) 



(3.4) 



c(A) d(X) 

where the m x m matrix functions b(X) and d(X) are given by 

6(A) = i9l(XI N - a)~%, d(X) = I m + i6* 2 {XI N -a)- 1 e 2 . (3.5) 
We first prove that 

b{X)d{X)- 1 = i6l{XI N - (3.6) 
Using ()1.3|) in preliminaries, from (j0.16|) and ()3.5|) we obtain 

d(X)- 1 = I m - i9* 2 (XI N - P)~%. (3.7) 
The equalities (jn~TH)l . (Q and (jSH|) yield 

b{X)d{X)- 1 = 

= i6l(XI N - a)- l 6 2 + i9l(XI N - aY\l3 - a)(XI N - py 1 9 2 . (3.8) 
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From (|3.8|) formula (J3.6)) follows. 

Let if be defined by (j(J.15|) . and thus by virtue of (|3.6j) we have 

tp{\) = b{X)d{\)- 1 . 



(3.9) 



By f!3.4|) . ()3.9|) and the representation ()0.10|) of the fundamental solution we 

get 



W n (X) 



I'm. 



X + i 
X 



W a>A (n, X) 





d(X)- 1 



(3.10) 



Notice also that (|U.14j) yields a more general formula than (J2.ll)) . namely: 

W a , A (n, X)*W a , A (n,X) = 

= I 2m - i(X - X)A* n (XI N - aT^-'i^N - a)- l X n . (3.11) 

As the second term in the right hand side of (|3.11|) is nonpositive, it follows 
from formula ()3.11)1 that 



W aA (n, XyW aA (n, A) < I 2m (A e C_). 



(3.12) 



Taking into account ()3.10|) and ()3.12|) we obtain (j(J.6p . i.e., tp is a Weyl 
function. 

It remains only to prove the uniqueness of the Weyl function. Let us first 
show that for some M > and all n > we have the inequality 



f*(i N - i(a*)-TKVN + i^Tf < Af/7, /el, 



(3.13) 



where 



L := span A0CT(Q) Im (Aijv - a) 
In view of (I0.13J) formula ()3.11|) yields 

, *\-l\nv-l 



9t(XI N -a 



In —i(a 



x 



x(I N + ia- l ) n (XI N - a)' 1 9 l < 



A- A 



/ m .(3.14) 



Now to get ()3.13|) from ()3.14|) we note that span^^Im (Aijv — a)~ 1 9i co- 
incides with the same span when A runs over an e-neighbourhood E of any 
Ao ^ &a, for any sufficiently small e > 0. 
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By (|3.11|) and (|3.13|) we can choose Mi > such that we have 



[I m 0)W a>A (n,\)*W a)A (n,X) 



1 m 




>- for all | A| > Mi. (3.15) 



Without loss of generality we may assume that Mi is large enough in order 
that Mi > \ \a\ \ and a(A) is invertible for |A| > Mi. Then, taking into account 
(HHU1) and (l3~T5J) . we obtain 



X 



xW n (A)*W n (A) 
for all A in the domain 



_ c(A)a(A) 



-i 



D = {A : |A| > Mi, 3 A < --}. 
In other words, for X E D we have 

oo 

^/W n (A)W n (A)/ = oo (feLi 

where 



>-(a(A)- 1 )*a(A)- 1 (3.16) 



1 



(3.17) 



n=0 



Ll '- lm c(A)a(A)- 1 

Suppose now that 99 and </3 are Weyl functions of (jU.4|) and that for some 
fixed Ao G D we have y5(A ) 7^ y?(A ). Put 



L? = Im 



+ Im 



Im 



According to the definition of the Weyl function we have 

00 

f*W n (\y W n (\) f <oc (/GL 2 ). 



(3.18) 



n=0 



As dimLi = m and dimL 2 > m, there is a non-zero vector / such that 
/ G (Li n L 2 ), which contradicts (l3~T7l) and (ETTSl . □ 

For the proof of Theorem IU.5I we shall use the following lemma which is 
of independent interest. 
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Lemma 3.2 A strictly proper rational m x m matrix function p admits a 
minimal realization of the form 



y{\)=i0l{\I N -(3)-% 



'2- 



(3.19) 



such that 



(3.20) 



Proof. We may assume that tp is given by the minimal realization (JU.17)) . 
First let us show that equation (j(J.18j) has a unique solution X > 0. 

The minimality of the realization (|U.17|) means that the pair {$*, 7} is 
observable and the pair {7, $2} is controllable. Notice that Im D Im 
and f*$$* = yields f*$ = 0, i.e., Im $ = Im #1?*. Hence the pair {7, 
^2^2} i s controllable too. Therefore the pair {"&2 n & 2 -> ^7*} i s observable. The 
pair {«7*, #1} is controllable and hence c-stabilizable, that is, there exists 
a matrix K such that 27* + d\K has all its eigenvalues in the open left 
half-plane. But then we can use Theorem 16.3.3 [21] (see also [20]) to show 
that the equation (|0.18J) has a unique non-negative solution X and that this 
solution X is positive definite. 

Next, let 0i, 2 , P be defined by (KTEty . From (KTTHJ) and (KTT91 it follows 
that ()3.20j) is satisfied. 

According to (|0.17|) and (j(J.19|) the function p is also given by the realiza- 
tion (|3.19|) . Moreover as the realization (j(J.17|) is minimal, the same is true 
for the realization (|3.19|) . □ 

Proof of Theorem IU.5L Let ip be a strictly proper rational m x m matrix 
function. Let 0i, 02, j3 be as in Lemma f3. 21 and put a = (3 + i0 2 02- Then the 
triple a, 0i, and 02 satisfies f|0.12|) . Furthermore, the pairs {(3, 82} and 
0i } are full range. Since a = (3 + Z0202, it is immediate that the pair {a, 62} 
is full range (see Section |TJ). By ()0.12j) we have 0* = a* + i9 2 9l = a - %Q\d\. 
Hence, as {/?*, 0i} is a full range pair, so the pair {a, 0i} is full range too. 
Therefore the triple a, 9\, and 02 is admissible. From Theorem ID .41 it follows 
now that the function {S n } determined by the admissible triple a, 9\ and 02 
is indeed a solution of the inverse problem. 

Let us prove now the uniqueness of the solution of the inverse problem. 
Suppose that there is system (j(J.4|) with another spin sequence {S n }, given by 
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the admissible triple a, 81, 62, and with the same Weyl function ip. According 
to Theorem IU.4I we have another realization for tp, namely 

l P (\) = i6* 1 (\Ik-(3)- 1 6 2 , p = a-ie 2 B*. (3.21) 

As the pairs {a, 81} and {a, 8 2 } are controllable, and (3 = a — i9 2 02, (3* = 
a — i9i9*, it follows that the pairs {(3, 82} and {f3*, 81} are also controllable. 
Thus the realization (|3.21|) is minimal and N = N. Moreover, there is (see 
Section ^) a state space similarity transforming the realization (JU.15)) into 
([3.21)1 . that is, there exists an invertible N x N matrix S such that 

(3 = S/3S- 1 , 8 2 = S8 2 , 9* x = 81S- 1 . (3.22) 

Recall that a, 81, 8 2 is an admissible triple, and therefore we have 

P-^* = l (e 1 d* l - 8 2 8l). (3.23) 

From (13~221 and (13~231) it follows that Z = S- 1 (S*)- 1 satisfies 

f3Z - Z(3* = i(Z8 1 8* 1 Z - 8 2 8* 2 ), Z > 0. (3.24) 

Completely similar to the uniqueness of X > in (|(J.18j) one obtains the 
uniqueness of the solution Z > of ()3.24|) . Compare now ()3.20|) and ()3.24|) 
to see that Z = In and thus S is unitary. In view of this we have 

a = SaS*, 8 2 = S8 2 , 0i = S6 l . 

The last transformation does not change the spin sequence, i.e., {S n } = {S n }. 
□ 

The case when i a (a) is important in the next section. We shall use 
the acronym FG to denote the class of spin sequences {S n } determined by 
the triples a, 8\, 82, with a an N x N non-singular matrix and 9\, 82 of size 
N x m, satisfying the identity ()0.12|) and the additional special condition 
i ^ a(a). Notice (see the beginning of the proof of Proposition l(J.2|) that 
(IU.12|) implies that a(a) C C+. The next lemma shows that without loss of 
generality we can also require that {a, 81} and {a, 8 2 } are full range pairs, 
i.e., FG C FG. 
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Lemma 3.3 Assume that the spin sequence {S n } belongs FG. Then it can 
be determined by a triple a, 9\, 62 such that a is non-singular, W.llty holds, 

1 ^ o~(a), and the pairs {a, 9{\ and {a, 62} are full range. 

Proof. Let N denote the minimal order of a (0, % ^ o~(a)) in the set of triples 
that satisfy (j0.12J) and determine the given spin sequence {S n }. Suppose the 
N x N matrix a and the N x m matrices 9\ and 62 form such a triple but 
the pair {a, 9 2 } is not full range. Put 

00 

L := span [^J lma k 9 2 , N : = dimL 

k=0 

and choose a unitary matrix q that maps L onto the L Q := Im[ijv 0]*. 
Then we have 

O2 ■= q02 = q , 

(3.25) 

where the A^o x N$ matrix 5 (0, i ^ cr(5)) and the No x m matrices 9\, 62 
form a triple which satisfies (|0.12|) and determines {S n }. To show this we 
need to make some preparations. As L is an invariant subspace of a, so 
L is an invariant subspace of a, and thus a has the block triangular form 
given in (|3.25j) . Moreover in view of the inclusion Im 62 C Lq, we have Im 

02 Q Lq, i.e., 62 has the block form given in ([3.25)1 . Taking into account 
that q is unitary 0, i cr(a) and that a, 61, 6 2 satisfy (|U.12|) and determine 
{S n }, we see that 0, i £ a(a) and that a, 9\, 9 2 satisfy ()U.12|) and determine 
{S n } too. So in view of ([3.25)1 we have 0, i $ a(a) and the triple a, 9\, 9 2 
satisfies (10.12)1 as well. Now use the matrices Q n defined in (j2.13[l to rewrite 
A*S~ 1 A n as a 2 x 2 block matrix with block of size m x m. This yields 

A n S n 1 A n = {(A n S n 1 A n )fcj}^ J - =1 



a := qaq 



a ct\2 

tt 2 2 



61 ■= q9 1 



9i 

K 
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with blocks 



(A;S; 1 A n ) 22 = 0;(J* + i(aTWtf-i(a*'~ J, ~'' 



^ l {I N + ia- l )- n {I N -ia- 1 ) n 6 2 , 

l\n( t ■( „ *\ — 1\ — n 1 /) 



(3.26) 



" 1 ) n ^2. 



Partition Qn+i - Qn into four blocks: Q n+ i -Q n = {Xkj}tj = i, where xn i 
an N x N block. In view of (KH3l . lETTfll) . and (|3~23J) we obtain 



is 



X21 = 0, X12 = 0, X22 = 0, Xn = 
x(I Nn -ia-Y9 2 e*(I No +i(a*) 



-2a~ l {I No + ioT' ) " 1 x 



-l\ni 



I N 



-l\-n-l 
-l\-n-l/ 



Denote by A n , S n , Q n , S n et cetera the matrices generated by the triple 5, 
6i, 6 2 . One can see that Xn — Qn+i — Qn- Taking into account Q = In 
we obtain that the matrices Q n are block diagonal: Q n = diag{Q n , In-n }- 
Now according to ()3.25j) . (|3.2fij) it follows that 



-A^S^ A n ~ A*S n 1 A n 



K*K 





(3.27) 



By (|U.9J) and (j3.27|) we get S n = S n , i.e., the triple 5, 6\, 9 2 determines the 
spin sequence {S n }. This contradicts the assumption that A" is minimal and 
therefore the pair {a, 9 2 } should be full range. 

In the same way we shall show that the pair {a, 9{\ is full range too. 
Indeed, suppose {a, 9{\ is not full range. Put now 

oo 

L := span [^J lma k 9x, N : = dimL 

fc=0 

and choose a unitary matrix q that maps Lq onto the Lq := Im [In 0]*. 
Then similar to the previous case we obtain 



qaq 



a CH12 
a 22 



q9i 



q9 2 



I7 2 

K 



(3.28) 
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where the N Q x N matrix a (0, i ^ cr(5)) and the N xm matrices 9i, 62 form 
a triple which satisfies (|(J.12|) . To show that the triple a, 9\, 62 determines 
S n we shall use the fact that a, 9 1 , 9 2 determines S n and rewrite A*E~ 1 A n 
clS cL 2 X 2 block matrix with blocks of size m x m as follows: 

h*rJ^n l ^n = { (A^£ n ^n) kj }k,j=l 

with 

(A*£~ A n ) 22 = 9lR~ 9 2 , 

(A;E; 1 A n ) 21 = 9*R-\l N -ia- 1 )- n ^N + ia- 1 ) n 9 1 , 

(A;E; 1 A n ) 12 = ^(/^-^-^(Jjv + ^a*)- 1 )-^ 1 ^, (3.29) 

(A^A^n = 9* 1 (I N -i(a*)- 1 ) n (lN + i(aT 1 )~ n Rn 1 x 

x (/jv - ia- l )- n (I N + icT 1 )^!. 

Here the matrices R n are given by (|2.12|) . Partition now R n+ \ — R n into four 
blocks: R n+ i — R n = {xkj}kj=v where xn is an No x block. In view of 
(EHoT) and (l3~2l?l) we obtain 

X21 = 0, X12 = 0, x 22 = 0, xu = ^-\l No - ia- l y n - 1 x 

x(/^ + ia-^MKlNo - i(5?*)- 1 ) n (/ i v + z^*)- 1 )-"-^^*)- 1 . 

Denote by A n , E n , i?„, 5 n et cetera the matrices generated by the triple 5, 
9i, 6> 2 . One can see that Xn — Rn+i — Rn- Taking into account R = In 
we obtain that the matrices R n are block diagonal: R n = diag{i? n , In~n }- 
Now according to (I3~2~K1) . it follows that 

By ()0.9|) and ()3.30|) we get S n = S n , i.e., the triple 5, Q\, 9 2 determines the 
spin sequence {S n }. So the pair {a, 9i} should also be full range. □ 

Finally, from the proof of Theorem IU.4I we have the following corollary. 

Corollary 3.4 Let the parameter matrices a (0, i o~{a)), E (0) > and 
A (0) satisfy the identity <\U. 7| ). Then the Weyl function if of the system 
determined by these matrices is given by the formula 





K*K 



(3.30) 
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This corollary is proved by transforming the matrices a, Eq(0) and A (0) 
into the equivalent set S 2 ckSq, In, S 2 A . 



4 Isotropic Heisenberg magnet 

Explicit solutions of the discrete integrable nonlinear equations form an in- 
teresting and actively studied domain (see references in [TJ [TTJ H21 1221 1213 
123 123 EH!)- To study the IHM model we insert an additional variable t in 
our notations: A n (£), £„(t), S n (t), W a< \(n, t, A), (p(t, A) and so on. Notice 
that the order N and the parameter matrix a do not depend on t, and the 
dependence on t of the other matrix functions is defined by the equations 

' /Ao = -2((a-iI N )- 1 A Q P + + (a + iI N )- 1 A P_), (4.1) 



dt 

P± = \(hm±J) 



dEo 
dt 



- ((a - iIn)- 1 ^) + (a + U N )- l E {t) + 

+E (t)(a* + UnY 1 + Z (t)(a* - iIn)" 1 + (4.2) 
+2(« 2 + In)" 1 (aA (t)JA (ty + A (t) JA (t)*a*) ((a*) 2 + I N )~ 



We assume that the parameter matrices a, £o(0) and Ao(0) satisfy the iden- 
tity 

a£ (t) " S (t)«* = zA (t)A (t)* (4.3) 

at t = 0. Then according to flOJ, (03) the identity ijOl) holds for all t. 
(This result can be obtained by differentiating both sides of ([4.3)1 .) 

Theorem 4.1 Assume the parameter matrices a (0, % (jL o~(a)), S (0) > 
and A (0) satisfy the identity 

aS (0) - E (0)a* = *A (0)A (0)*. 

Define S (t) and A (t) by equations fl^.ip and Q^.ifl . T/ien X (£) > on 
some interval —e < t < e, and the sequence {S n (t)} given by (j£A,9j) and (jfl.&j) 
belongs to FG for each t from this interval. Moreover, {S n (t)} (—e < t < e) 
satisfies the IHM equations (|Q3ij) . (HOg). 
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Proof. As a does not depend on t and (|4.3|) is true, it is immediate that 
{S n (t)} c FG. 

Similar to the cases treated in jSU E2] we shall successively obtain the 
derivatives 4j.A n , 4£ n ) ^(A^E" 1 ), and ^W a ^(n, t, A), and use the expres- 
sions for these derivatives to derive the zero curvature equation ()0.20|) , which 
is equivalent to (jTHH|) . (JTHHjl : see [12]. In view of (jtTHIjl and flUH) we have 

A n (t) = [(I N + ia- 1 ) n e- 2t{a - ilN) ~ 1 e l (I N - ia- 1 ) n e- 2tia+ilN) ~ 1 6 2 ] , (4.4) 

where 9 P := (A (0)) p (p = 1,2). Hence it follows that 

^ = -2 ((a - U N y l K n P + + (a + il^KP-)- (4.5) 
Now we shall show by induction that 

^ = -((a-iI N y 1 Z n (t) + (a + iI N y 1 X n (t) + 

+S n (t)(a* + z/^)- 1 + E n (t)(a* - Un)" 1 + (4.6) 
+2(a 2 + / w )" 1 (aA n (t)JA n (t)* + A n (t) JA n (t)*a*) ((a*) 2 + In)- 1 ). 



By ()4.2|) formula ()4.6|) is true for n = 0. Suppose it is true for n = r. Then, 
taking into account the second relation in (j().8j) we obtain 

dSr+1 - ^(a-2/ J v)- 1 S r+1 (t)-(a-2/ JV )- 1 a- 1 A r (t)JA r (t)*(a*)-N 



+ {a + Un)- 1 ^^) - {a + U N y l a- l K(t) JA^i)*^*)- 1 + 
+£ r+ i(*)(a* + iI N )~ l - a- 1 K r {t)JK{ty{a*)-\a f + iI N y l + 
+£ r+ i(*)(a* - z/^v)' 1 - a- 1 A r {t)JA r {ty{a*y 1 {a* - iIn)' 1 + 
+2(a 2 + I N y\(a + a- l )A r (t)JA r (ty + A r (t)JA r (ty x (4.7) 

x(a* + (a*)- 1 )) ((a*) 2 + J^)" 1 ) + | {a~ l A r {t) J A r (t)* (a*y l ) . 
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(4.8) 



+(q - il^a^A^t) JA r (*)*(a*) _1 + (a + iI N )~ 1 a~ 1 A r (t)J x 

xA r (£)*(a*) -1 + JA r (£)*(a*)~V* + U N y l + a" 1 x 

xA r (t) JA r {t)*{a*)-\a* - U N y l = -(a - il N )~ l a~ l A r (t) x 
x\ r {ty(a*y l + (a + z/ JV )- 1 a- 1 A r (t)A r (t)*(a*)- 1 - cr^t) x 



Notice that (a — Un) 1 — (a + z/jv) 1 = 2z(a 2 + Jjv) x . Therefore we can 
rewrite ()4.8|) as 



d(t) = 2i(a- 1 A r (t)A r (ty(a*y 1 ((u*) 2 + I N y 1 - 



C 2 (t) : = 2(a 2 + In)" 1 ((a + a _1 )A r (t) JA r (t)* + A r (t)JA r (t)* x 
x(a* + (a*)" 1 )) ((a*) 2 + J^)- 1 = 2 JA r (t)* x 
x((a*) 2 + I^- 1 + (a 2 + JA^V*)" 1 ). (4.10) 



Using the first relation in (j(J.8|) and equalities (|4.9|) . (|4.1U|) we get 



C 2 (t) - C x {t) = 2{a- l A r {t)JA r+l {ty{{a*f + I N y l + 

+ (a 2 + I^Ar+^JArityiaT 1 ) = 2(a 2 + In)' 1 x 
x (aA r+1 (t) JA r+1 (t)* + A r+1 (t) JA r+1 (t)*a*) ((a*) 2 + I N y\ (4.11) 



From (|4.7jl and (|4.11|) it follows that (jUBj) is valid for n = r + 1 and so for 



Taking into account (J4.5|) and (J4.6j) we can obtain the equation 



-(A,(t)*E„(t)- 1 ) = //+(t)A n (i)*E n (i)- 1 ( Q; -a iV )- 1 + 



(a 2 + /at) 1 a 



(4.9) 



Notice also that 



all n > 0. 



+#-(t)A n (t)*E n (£)-V + 2/ JV ) 



-l 



(4.12) 



where 



#+(t) = 2W aA (n, t, i)P+W aA (n, t, -«)*, 



(4.13) 
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H~(t) = 2W a>A (n,t,-i)P-W a>A (n,t,i)*. 



(4.14) 



Indeed, by ()4.5|) we have 



-( k (tyj: n (t)- 1 ) = -(2P + A n (ty(a* + ti N )- 1 + 

+2P_A n (t)>* - iIm)- 1 + E n (t)- 1 ^(t))E n (t) 



(4.15) 



Identity (UTT^ yields 



+z((a*) ± 2/ JV )' 1 S„(t)- 1 A n (t)A„(t)*S„(t)- 1 ( a ± z/^)- 1 . (4.16) 



By using f)4.6|) . ()4.16|) and (|4.17|1 . and after some calculations, we rewrite 
(j4.15j) in the form (|4.12[) . where 



= J 2m + J-^A„(t)*S„(t)- 1 (a-^)- 1 A n (t)J + 

+zJA n (t)*(a* - i/ JV )- 1 E n (i)- 1 A n (i) + A„(t)*E n (t)- 1 x 
x(o - z/ iV )~ 1 A n (t)JA n (t)*(a* - iJjv)- 1 ^,,^)- 1 ^^), (4.18) 

#~(t) = I 2m - J + ik n (ty^ n (t)-\a + iI N )- l k n (t)J - 

-iJA n (t)*(a* + z/ JV )- 1 E n (t)- 1 A n (t) - A n (t)*E„(t)- 1 x 
x(a + z/ JV )- 1 A n (t)/AnW*(a* + ^iv)~ 1 S n (t)- 1 A n (t). (4.19) 

From (jnXTJ) and (gHHD it follows that 



and so, taking into account (|2.11|) . we derive (|4.13|) . Equality (|4.14j) follows 
from ()4.19|) in a similar way. 



Finally notice that 



2(a 2 + I N )~ l a ={a- iI N )~ l + (a + U N ) 



-i 



(4.17) 



H+(t) = I 2m + W aA (n, t, i)JW aA (n, t, -i)* 
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Recall now that m = 2 and (l2~TTf holds. Then according to fl2T7p , (|2~2Tf 
and (|4.13|) we get 

H+(t) = c+(t)(I 2 + S n (t))(I 2 + S n ^(t)), (4.20) 
and according to (|2~T%1) . (I2~2"2]) and (ETT41 we get 

H-(t) = c-(t){I 2 - S n (t)){I 2 - S n ^(t)), (4.21) 

where c^(t) are scalar functions. In view of (j4.13j) and (j4.14|) we obtain also 
Tr H^, where Tr denotes the trace. Indeed, 

Tr#±(t)=2. (4.22) 

By Remark 12.41 and formula (|U.23|) we derive 

Tr (I 2 ±S n (t))(I 2 ±S n ^(t) ) = Tr (I 2 + S n (t)S n ^(t)) = 

= 2(1+ 5 n _!(i) • ~S n (t)). (4.23) 

From dO0|)-(jO3J) it follows that 1 + SVi(0 • S\^) ^ 0. Now formulas 
(HT22J), flOgD and (l4"22j) yield that 

-&^(t) = 2 = '&fl*(0. (4.24) 

Taking into account (KT22j) . (jj^Uj) and (ET2TD we see that V± = and 
so ()4.24|) yields equalities = if*. Hence we have 

^(A^)*^)- 1 ) = V+(t)A n {t)*Z n {t)- 1 {a-iI N )- 1 + 

+V-(t)A n (t)*Z n (t)- 1 (a + iI N )~ l . (4.25) 

From V± = H±, pTT]) and definitions iftTHijl and (jglj) we also get: 

^W^Afn, t, ±«) = 2W aA {n, t, ±i)P ± . (4.26) 

Let us differentiate now W a ^. For this purpose notice that 

I _i_ • j \-\( XT \-i (a-A/jv)~ 1 -(Q±2/j V )~ 1 

(a±z/jv) (a-A/jv) = r— — • (4.27) 

A ± % 
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Using (jOl) . $£2B and f *^l ) w e derive 

d nr t * \\ t/+/.\ W a A(n, t, A) — W a \(n, t, i) 

— W aA (n, t, A) = V+(t) ; — + V n (t) x 

W«,a(M,A) - W ai A(M,-«) n^A^U) - Wa,A(«,*,0 D 

x TT^ z \ : + 

A ■ / A — 2 

2 W a , A (n,t, A) - Wa, A (w,t,-z) p 2g 

X + i " [ ' } 

In view of (14. 26(1 we can rewrite (|4.28jl as 

^-W a , A (n, t, A) = F n (t, A) W aA (n, t, A) - W aA (n, t, X)F(t, A), (4.29) 
at 

where F n is given by the second relation in (j0.21|) and 

F = 2((\-i)- 1 P + + (A + i)- x P_). 

Thus in view of Theorem IU.1I and formula ()4.29|) the non-degenerate matrix 
functions {W n } given by 



W n (t,\) = \- n W aA (n,t,\) 



(A - iye^-i)- 1 o 
(A + iYe^+i)- 1 



satisfy equations (|0.25j) . i.e., the compatibility condition ()0.20|) is valid. As 
equation (j(J.20|) is equivalent to ()0.24|) . the theorem is proved. □ 

Theorem 14. II together with Corollary 13.41 yields the following result. 

Corollary 4.2 Under the conditions of Theorem \4-l\ the evolution of the 
Weyl function tp of the system W n +i(t, A) = G n (t, X)W n (t, A) is given by the 
formula 

<p(t, A) = 2^e- 2 * (Q * + ^ rl S (t)- 1 (AJ A . - (3(t))- 1 e- 2t ^ a+lI ^' 1 9 2 , (4.30) 

pit) = a - ie- 2t{a+ilNrl e 2 e* 2 e~ 2t ^- ilN ^ 1 j: (t)- 1 . 

As an illustration let us consider a simple example. 

Example 4.3 Put m — n — 1 and a = ih (h > 0, h ^ 1), and choose 
scalars 9\, 9 2 such that \9i\ 2 + \9 2 \ 2 = Ih. Then a, 9\, 9 2 form an admissible 
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triple and a, S (0) = 1, A (0) = [Q\ 9 2 } satisfy the conditions of Theorem 
\4-l\ and Corollary \4-&\ Therefore by ( |^.^| ) we have 

A n {t) = h- n [{h + l) n 6 1 e X p{^} (^-lr^expj^-}]. (4.31) 

From <\U.14\) and ^~JJl it follows that 



E„(f) = c n (h) := (h + l) 2 "^ 2 + (h- l) 2n \9 2 \ 2 . (4.32) 



According to (pTffi . ( P~37] ) 7 and we get 

shy^Wh 2 -l) 2n 

c n (h)c n+1 (h) 



now 



(5 B (t)) U = 1 ' 1 * K , u ' , (Sn(t))22 = -(Sn(t))ll, 



,„ . ,„ . Ah9 x 9 2 r 4z£ 1 

(*(t)) u = (*(t))i - ^CT exp {r^} x 

x(/i 2 - l) n ((/i + l) 2n+1 |^| 2 -(h- l) 2n+1 |6» 2 | 2 ). 
Finally Corollary \4.ty yields: 

Ait 1 ^6> 2 



p(t,A) = exp j— 



h 2 ) \ + i{\e 2 \ 2 - h) 
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